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The Mott and the Anderson routes to locahzation have long been recognized as the two basic 
processes that can drive the metal-insulator transition (MIT). Theories separately describing each 
of these mechanisms were discussed long ago, but an accepted approach that can include both 
has remained elusive. The lack of any obvious static symmetry distinguishing the metal from the 
insulator poses another fundamental problem, since an appropriate static order parameter cannot 
be easily found. More recent work, however, has revisited the original arguments of Anderson 
and Mott, which stressed that the key diference between the metal end the insulator lies in the 
dynamics of the electron. This physical picture has suggested that the "typical" (geometrically 
averaged) escape rate = exp(lnrj^c) from a given lattice site should be regarded as the proper 
dynamical order parameter for the MIT, one that can naturally describe both the Anderson and the 
Mott mechanism for localization. This article provides an overview of the recent results obtained 
from the corresponding Typical- Medium Theory, which provided new insight into the the two-fluid 
character of the Mott- Anderson transition. 

I. FROM METAL TO INSULATOR: A NEW PERSPECTIVE 

Metal or insulator - and why? To answer this simple question has been the goal and the driving force for 
much of the physical science as we know it today. Going back to Newton's not-so-successful exercises in Alchemy, 
the scientist had tried to understand what controls the flow of electricity in metals and what prevents it in insulators 
[1]. To understand it and to control it - achieving this could prove more useful and lucrative than converting lead 
into gold. Indeed, the last few decades have witnessed some most amazing and unexpected advances in material 
science and technology. And this ability - its intellectual underpinning - is what was indispensable in designing and 
fabricating the iPhone, the X-Box, and the MRI diagnostic tool. Today's kids have grown up in a different world 
than had their parents - all because we have learned a few basic ideas and principles of electron dynamics. 

In almost every instance, these advances are based on materials that find themselves somewhere between metals 
and insulators. Material properties are easy to tune in this regime, where several possible ground states compete [S]. 
Here most physical quantities display unusual behavior and prove difficult to interpret using conventional ideas 
and approaches. Over the last few decades, scores of theoretical scenarios and physical pictures have been proposed, 
most of which will undoubtedly end up in back drawers of history. Last couple of years, however, have seen a veritable 
avalanche of new experimental results, which provide compelling clues as to what the theorists should not overlook: 
the significant effects of spatial inhomogeneities in the midst of strongly correlated phases. 




Figure 1: Spectacular advances in scanning tunneling microscopy (STM) have revealed that many "bad metals" or barely- 
doped insulators are surprisingly inhomogeneous on the nano-scale. Understanding and controlling these materials will not 
be possible without coming to grips with the origin, the stability, and the statistics of such mesoscopic granularity. (A) 
"Tunneling asymmetry" imaging 2 provides evidence for the emergence of a low temperature "electronic cluster glass" within 
the superconducting phase of Cai,^g,NaQ,i2Cu02Cl2', (B) Fourier-transform STM "15] reveals nano-scale Fermi surface variations 
in Bi2Sr2CuOe-x; (C) Differential conductance maps \4l showing spatial variations of the local pseudogaps in the normal phase 

(T > Te) of Bi2Sr2CaCu20s+5. 
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To understand many, if not most exotic new materials, one has to tackle the difficult problem of understanding 
the metal-insulator quantum phase transition, as driven by the combined effects of strong electronic correlations 
and disorder. Traditional approaches to the problem, which emerged in the early 1980s, have focused on examining 
the perturbative effects of disorder within the Fermi liquid framework. Despite their mathematical elegance, these 
theories, unfortunately, prove ill-suited to describe several key physical processes, such as tendency to local magnetic 
moment formation and the approach to the Mott insulating state. In addition, such weak-coupling theories cannot 
easily describe strongly inhomogeneous phases, with behavior often dominated by broad distributions and rare disorder 
configurations. 

This new insight, which is largely driven by experimental advances, seems to suggest that an alternative theoretical 
picture may provide a better starting point. In this article we describe recent advances based on a new theoretical 
method, which offers a complementary perspective to the conventional weak-coupling theories. By revisiting the 
original ides of Anderson and Mott, it examines the typical escape rate from a given site as the fundamental dynamical 
order parameter to distinguish between a metal and an insulator. This article describes the corresponding Typical- 
Medium Theory (TMT) and discusses some of its recent results and potential applications. We fist discuss discuss, 
in some detail, several experimental and theoretical clues suggesting that a new theoretical paradigm is needed. The 
formulation of TMT for Anderson localization of noninteracting electrons is then discussed, with emphasis on available 
analytical results. Finally, we review recent progress in applying TMT to the Mott- Anderson transition for disordered 
Hubbard models, and discuss resulting the two-fluid behavior at the critical point. 

II. THEORETICAL CHALLENGES: BEYOND CINDERELLA'S SLIPPER? 

The existence of a sharp metal- insulator transition at T = has been appreciated for many years [T] . Experiments 
on many systems indeed have demonstrated that a well defined critical carrier concentration can easily be identified. 
On the theoretical side, ambiguities on how to describe or even think of the metal-insulator transition have made it 
difficult to directly address the nature of the critical region. In practice, one often employs the theoretical tools that 
are available, even if possibly inappropriate. Even worse, one often focuses on those systems and phenomena that fit 
an available theoretical mold, ignoring and brushing aside precisely those features that seem difficult to understand. 
This "Cinderella's slipper" approach is exactly what one should not do; unfortunately it happens all too often. A cure 
is, of course, given by soberly confronting the experimental reality: what seems paradoxical at first sight often proves 
to be the first clue to the solution. 



A. Traditional approaches to disordered interacting electrons 

Most studies carried out over the last thirty years have focused on the limit of weak disorder [7] , where consider- 
able progress has been achieved. Here, for non-interacting electrons the conductance was found to acquire singular 
(diverging) corrections in one and two dimensions, an effect known as "weak localization" [7l . According to these pre- 
dictions, for d <2 the conductivity would monotonically decrease as the temperature is lowered, and would ultimately 
lead to an insulating state at T = 0. Interestingly, similar behavior was known in Heisenberg magnets [5Mll|. where 
it resulted from d = 2 being the lower critical dimension for the problem. This analogy with conventional critical 
phenomena was first emphasized by the "gang of four" [8] , as well as Wegner |9l [10] , who proposed an approach to 
the metal-insulator transition based on expanding around two dimensions. For this purpose, an effective low energy 
description was constructed [HI [101 [H] i which selects those processes that give the leading corrections at weak disorder 
in and near two dimensions. This "non- linear sigma model" formulation [3 \W\ [T^] was subsequently generalized 
to interacting electrons by Finkelshtein [13_, and studied using renormalization group methods in 2 -I- e dimensions 
[13-15 . In recent years, the non-linear sigma model of disordered interacting electrons has been extensively studied 
by several authors |TFl [T7| . 

While the sigma model approach presented considerable formal complexity, its physical content proved - in fact 
- to be remarkably simple. As emphasized by Castellani, Kotliar and Lee (TS], one can think of the sigma model 
of disordered interacting electrons as a low energy Fermi liquid description of the system. Here, the low energy 
excitations are viewed as a gas of diluted quasi-particles that, at least for weak disorder, can be described by a small 
number of Fermi liquid parameters such as the diffusion constant, the effective mass, and the interaction amplitudes. 
In this approach, one investigates the evolution of these Fermi liquid parameters as weak disorder is introduced. The 
metal-insulator transition is then identified by the instability of this Fermi liquid description, which in d = 2 -|- e 
dimensions happens at weak disorder, where controlled perturbative calculations can be carried out. 
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Remarkably, by focusing on such a stability analysis of the metallic state, one can develop a theory for the transition 
which does not require an order parameter description, in contrast to the standard approaches to critical phenomena 
This is a crucial advantage of the sigma model approach, precisely because of the ambiguities in defining 
an appropriate order parameter. We should stress, however, that by construction, the sigma model focuses on 
those physical processes that dominate the perturbative, weak disorder regime. In real systems, the metal-insulator 
transition is found at strong disorder, where a completely different set of processes could be at play. 

B. Anderson's legacy: strong disorder fluctuations 

From a more general point of view, one may wonder how pronounced are the effects of disorder on quantum phase 
transitions. Impurities and defects are present in every sample, but their full impact has long remained ill-understood. 
In early work, Griffiths discovered [19} that rare events due to certain types of disorder can produce nonanalytic 
corrections in thermodynamic response. Still, for classical models and thermal phase transitions he considered, these 
effects are so weak to remain unobservably small [2Ij. The critical behavior then remains essentially unmodified. 

More recent efforts turned to quantum (T = 0) phase transitions [22) , where the rare disorder configurations prove 
much more important. In some systems they gives rise to "Quantum Griffiths Phases" (QGP) |5J [53], associated 
with the "Infinite Randomness Fixed Point " (IRFP) phenomenology!^. Here, disorder effects produce singular 
thermodynamic response not only at the critical point, but over an entire region in its vicinity. In other cases, related 
disorder effects are predicted |251 E5] to result in "rounding" of the critical point, or to produce intermediate "cluster 
glass" phases |27l [25] masking the critical point. Physically, QGP-IRFP behavior means (5] [53] that very close to the 
critical point, the system looks increasingly inhomogeneous even in static response. 

But how robust and generic may such pronounced sensitivity to disorder be in real systems? Does it apply only 
to (magnetic and/or charge) ordering transitions, or is it relevant also for the metal-insulator transitions (MITs)? A 
conclusive answer to these questions begs the ability to locally visualise the system on the nano-scale. Remarkably, 
very recent STM images provide striking evidence of dramatic spatial inhomogeneities in surprisingly many systems. 
While much more careful experimental and theoretical work is called for, these new insights makes it abundantly 
clear that strong disorder effects - as first emphasized by early seminal work of Anderson [35] - simply cannot be 
disregarded. 




0.4 1.0 1.4 



Figure 2: Theory predictions [29] for an "Electronic Griffiths Phase" [30; in a moderately disordered normal metal near a Mott 
meal-insulator transition. Local density of states (LDOS) spectra look dramatically "smoother" near the Fermi energy (bottom) 
than away from it (top). This contrast is more pronounced close to the Mott transition (right), than outside the critical region 
(left). Very similar behavior was recently observed by STM imaging of the superconducting phase of doped cuprates ,2,, but our 
results strongly suggest that such energy-resolved "disorder healing" |29H31] is a much more general property of Mott systems. 
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C. The curse of Mottness: the not-so-Fermi liquids 

One more issue poses a major theoretical challenge. According to Landau's Fermi liquid theory, any low temperature 
metal behaves in a way very similar to a gas of weakly interacting fermions. In strongly correlated systems, closer to 
the Mott insulating state, this behavior is typically observed only below a modest crossover temperature T* <C Tp- 
Adding disorder typically reduces T* even further, and much of the experimentally relevant temperature range simply 
does not conform to Landau's predictions. Theoretically, this situation poses a serious problem, since the excitations 
in this regime no longer assume the character of diluted quasiparticles. Here perturbative corrections to Fermi liquid 
theory simply do not work [B], and a conceptually new approach is needed. 
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Figure 3: Finite temperature metal-insulator crossover in transport close to a disordered Mott transition. Very high values of 
resistivity, strongly exceeding the "Mott limit" [T] |33] are observed above the crossover temperature T* . Remarkably, increasing 
disorder W reduces the resistivity maximum, rendering the system effectively more metallic. This behavior, which is clearly seen 
in our DMFT modelling 1341 (left panel), has very recently been also observed in experiments pSj on organic charge-transfer 
salts (right panel), where disorder is systematically introduced by X-ray irradiation. 

A new theoretical paradigm, which works best precisely in the incoherent metallic regime, has been provided by the 
recently developed Dynamical Mean-Field Theory (DMFT) methods [36 . Unfortunately, in its original formulation, 
which is strictly exact in the limit of infinite dimensions, DMFT is not able to capture Anderson localization effects. 
Over the last twelve years, this nonperturbative approach has been further extended jHl 1571 - lii] to incorporate the 
interplay between the two fundamental mechanisms for electron localization: the Mott (interaction-driven) [IJ and 
the Anderson (disorder-driven) [321 US] route to arrest the electronic motion. In addition, the DMFT formulation can 
be very naturally extended to also describe strongly inhomogeneous and glassy phases of electrons ^46n55j , and even 
capture some aspects of the Quantum Griffiths Phase physics found at strong disorder jBl HH HHJ ISH HOI HH . 
In the following, we first discuss the DMFT method as a general order-parameter theory for the metal-insulator 
transition, and the explain how it needs to be modified to capture Anderson localization effects. 

III. ORDER-PARAMETER APPROACH TO INTERACTION-LOCALISATION 

A. Need for an order-parameter theory: experimental clues 

In conventional critical phenomena, simple mean-field approaches such as the Bragg- Williams theory of magnetism, 
or the Van der Waals theory for liquids and gases work remarkably well - everywhere except in a very arrow critical 
region. Here, effects of long wavelength fluctuations emerge that modify the critical behavior, and its description 
requires more sophisticated theoretical tools, based on renormalization group (RG) methods. A basic question then 
emerges when looking at experiments: is a given phenomenon a manifestation of some underlying mean-field (local) 
physics, or is it dominated by long-distance correlations, thus requiring an RG description? For conventional criticality 
the answer is well know, but how about metal-insulator transitions? Here the experimental evidence is much more 
limited, but we would like to emphasize a few well-documented examples which stand out. 
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1. Doped semiconductors 



Doped semiconductors such as Si:P [SJ are the most carefully studied examples of the MIT critical behavior. Here 
the density-dependent conductivity extrapolated to T = shows sharp critical behavior [G^ of the form a ^ {n — Uc)^, 
where the critical exponent « 1/2 for uncompensated samples (half-filled impurity band), while dramatically 
different w 1 is found for heavily compensated samples of Si:P,B, or in presence of strong magnetic fields. Most 
remarkably, the dramatically differences between these cases is seen over an extremely broad concentration range, 
roughly up to several times the critical density. Such robust behavior, together with simple apparent values for the 
critical exponents, seems reminiscent of standard mean-field behavior in ordinary criticality. 




Figure 4: Critical behavior of the conductivity for uncompensated Si:P and compensated Si:P,B |65| . The conductivity exponent 
~ 1/2 in absence of compensation, while ~ 1 in its presence. Clearly distinct behavior is observed in a surprisingly broad 
range of densities, suggesting mean-field scaling. Since compensation essentially corresponds to carrier doping away from a 
half-filled impurity band [64j , it has been suggested [7] that the difference between the two cases may reflect the role of strong 
correlations. 



2. 2D-MIT 



Signatures of a remarkably sharp metal-insulator transition has also been observed |66fl68| in several examples of 
two-dimensional electron gases (2DEG) such as silicon MOSFETs. While some controversy regarding the nature or 
even the driving force for this transition remains a subject of intense debate, several experimental features seem robust 
properties common to most studied samples and materials. In particular, various experimental groups have demon- 
strated [661 157] striking scaling of the resistivity curves in the critical region, which seems to display [BS] remarkable 
mirror symmetry ("duality") |70l over a surprisingly broad interval of parameters. In addition, the characteristic 
behavior extends to remarkably high temperatures, which are typically comparable the Fermi temperature [68]. One 
generally does not expect a Fermi liquid picture of diluted quasiparticles to apply at such "high energies", or any 
correlation length associated with quantum criticality to remain long. 

These experiments taken together provide strong hints that an appropriate mean-field description is what is needed. 
It should provide the equivalent of the Van der Waals equation of state, for disordered interacting electrons. Such a 
theory has long been elusive, primarily due to a lack of a simple order-parameter formulation for this problem. Very 
recently, an alternative approach to the problem of disordered interacting electrons has been formulated, based on 
dynamical mean-field (DMFT) methods PB]. This formulation is largely complementary to the scaling approach, and 
has already resulting in several striking predictions. In the following, we briefly describe this method, and summarize 
the main results that have been obtained so far. 
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Figure 5: The resistivity curves (left panel) for a two-dimensional electron system in silicon [66, show a dramatic metal-insulator 
crossover as the density is reduced below ric ~ W^^ cm~^ . Note that the system has "made up its mind" whether to be a metal 
or an insulator even at surprisingly high temperatures T ~ Tf ~ lOK. The right panel displays the scaling behavior which 
seems to hold over a comparable temperature range. The remarkable "mirror symmetry" I69j of the scaling curves seems to hold 
over more then an order of magnitude for the resistivity ratio. This surprising behavior has been interpreted [70j as evidence 
that the transition region is dominated by strong coupling effects characterizing the insulating phase. 



B. The DMFT physical picture 

The main idea of the DMFT approach is in principle very close to the original Bragg- Williams (BW) mean-field 
theories of magnetism |71] . It focuses on a single lattice site, but replaces |36| its environment by a self-consistently 
determined "effective medium", as sho"wn in Fig. 1.3. 




Figure 6: In dynamical mean-field theory, the environment of a given site is represented by an effective medium, represented by 
its "cavity spectral function" Ai(cj). In a disordered system, Ai(aj) for different sites can be very different, reflecting Anderson 
localization effects. 

In contrast to the BW theory, the environment cannot be represented by a static external field, but instead must 
contain the information about the dynamics of an electron moving in or out of the given site. Such a description 
can be made precise by formally integrating out all the degrees of freedom on other lattice sites. In presence of 
electron-electron interactions, the resulting local effective action has an arbitrarily complicated form. Within DMFT, 
the situation simplifies, and all the information about the environment is contained in the local single particle spectral 
function Aj(w). The calculation then reduces to solving an appropriate quantum impurity problem supplemented by 
an additional self-consistency condition that determines this "cavity function" Ai(aj). 
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The precise form of the DMFT equations depends on the particular model of interacting electrons and/or the form 
of disorder, but most applications ^36j to this date have focused on Hubbard and Anderson lattice models. The 
approach has been very successful in examining the vicinity of the Mott transition in clean systems, in which it 
has met spectacular successes in elucidating various properties of several transition metal oxides |40j , heavy fermion 
systems, and even Kondo insulators |72J. 

C. DMFT as an order-parameter theory for the MIT 

The central quantity in the DMFT approach is the local "cavity" spectral function Ai(a;). From the physical point 
of view, this object essentially represents the available electronic states to which an electron can "jump" on its way 
out of a given lattice site. As such, it provides a natural order parameter description for the MIT. Of course, its form 
can be substantially modified by either the electron-electron interactions or disorder, refiecting the corresponding 
modifications of the electron dynamics. According to Fermi's golden rule, the transition rate to a neighboring site 
is proportional to the density of final states - leading to insulating behavior whenever Ai{uj) has a gap at the Fermi 
energy. In the case of a Mott transition in the absence of disorder, such a gap is a direct consequence of the strong 
on-site Coulomb repulsion, and is the same for every lattice site. 
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Figure 7: The local cavity spectral function Ai{ij) as the order parameter for the MIT. In a metal (a) there are available 
electronic states near the Fermi level (dashed line) to which an electron from a given site can delocalize. Both for a Mott 
insulator (b) and the Anderson insulator (c) the Fermi level is in the gap, and the electron cannot leave the site. Note that 
the averaged spectral function (dotted line in (c)) has no gap for the Anderson insulator, and thus cannot serve as an order 
parameter. 

The situation is more subtle in the case of disorder-induced localization, as first noted in the pioneering work of 
Anderson Here, the average value of Ai{uj) has no gap and thus cannot serve as an order parameter. However, 
as Anderson noted a long time ago, "...no real atom is an average atom..." |45| . Indeed, in an Anderson insulator, 
the environment "seen" by an electron on a given site can be very different from its average value. In this case, the 
typical "cavity" spectral function Ai^ui) consists of several delta- function (sharp) peaks, reflecting the existence of 
localized (bound) electronic states, as shown in Fig. 1.4(c). Thus a typical site is embedded in an environment that 
has a gap at the Fermi energy - resulting in insulating behavior. We emphasize that the location and width of these 
gaps strongly vary from site to site. These strong fluctuations of the local spectral functions persist on the metallic 
side of the transition, where the typical spectral density Atyp = exp < In(Ai) > can be much smaller than its average 
value. Clearly, a full distribution function is needed to characterize the system. The situation is similar as in other 
disordered systems, such as spin glasses [73 . Instead of simple averages, here the entire distribution function plays a 
role of an order parameter, and undergoes a qualitative change at the phase transition. 

The DMFT formulation thus naturally introduces self-consistently deflned order parameters that can be utilized to 
characterize the qualitative differences between various phases. In contrast to clean systems, these order parameters 
have a character of distribution functions, which change their qualitative form as we go from the normal metal to the 
non-Fermi liquid metal, to the insulator. 
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IV. TYPICAL MEDIUM THEORY FOR ANDERSON LOCALIZATION 



In the following, we demonstrate how an appropriate local order parameter can be defined and self-consistently 
calculated, producing a mean-field like description of Anderson localization. This formulation is not restricted to 
either low temperatures or to Fermi liquid regimes, and in addition can be straightforwardly combined with well- 
known dynamical mean-field theories (DMFT) [551 [551 IT^HTS] of strong correlation. In this way, our approach 
which we call the Typical Medium Theory (TMT), opens an avenue for addressing questions difficult to tackle by any 
alternative formulation, but which are of crucial importance for many physical systems of current interest. 

Our starting point is motivated by the original formulation of Anderson |32| . which adopts a local point of view, 
and investigates the possibility for an electron to delocalize from a given site at large disorder. This is most easily 
accomplished by concentrating on the (unaveraged) local density of electronic states (LDOS) 

=^5(c^-c^„)|^„(i)|2. (1) 

n 

In contrast to the global (averaged) density of states (ADOS) which is not critical at the Anderson transition, the 
LDOS undergoes a qualitative change upon localization, as first noted by Anderson [35]. This follows from the fact 
that LDOS directly measures the local amplitude of the electronic wavefunction. As the electrons localize, the local 
spectrum turns from a continuous to an essentially discrete one [35], but the typical value of the LDOS vanishes. 
Just on the metallic side, but very close to the transition, these delta-function peaks turn into long-lived resonance 
states and thus acquire a finite escape rate from a given site. According to to Fermi's golden rule, this escape rate 
can be estimated [35] as ^ t^p, where t is the inter-site hopping element, and p is the density of local states of 
the immediate neighborhood of a given site. 

The typical escape rate is thus determined by the typical local density of states (TDOS), so that the TDOS directly 
determines the conductivity of the electrons. This simple argument strongly suggests that the TDOS should be 
recognized as an appropriate order parameter at the Anderson transition. Because the relevant distribution function 
for the LDOS becomes increasingly broad as the transition is approached, the desired typical value is well represented 
by the geometric average ptyp = exp{< Inp >}. Interestingly, recent scaling analyses [771 178] of the multi-fractal 
behavior of electronic wavefunctions near the Anderson transition has independently arrived at the same conclusion, 
identifying the TDOS as defined by the geometric average as the fundamental order parameter. 



A. Self-consistency conditions 



To formulate a self-consistent theory for our order parameter, we follow the "cavity method," a general strategy that 
we borrow from the DMFT [36]. In this approach, a given site is viewed as being embedded in an effective medium 
characterized by a local self energy function S(i^). For simplicity, we concentrate on a single band tight binding 
model of noninteracting electrons with random site energies Si with a given distribution P{ei). The corresponding 
local Green's function then takes the form 



G{Lo,e,) = [Lo-e,- ^{uj)]-\ (2) 

A(a;) = Ao(a;-E(a;)) = A'+iA", (3) 

A„(w) =c^-l/Go(cj), (4) 

G„(c.)= r^d.'^ (5) 

is the Hilbert transform of the bare density of states po{oj) which specifies the band structure. 

Given the effective medium specified by a self-energy S(a;), we are now in the position to evaluate the order 
parameter, which we choose to be the TDOS as given by 

Ptyp(w) = cxp / dsi P{ei) lnp{uj,ei) \ , (6) 



Here, the "cavity function" is given by 



and 



where the lattice Green's function 
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^lm.G{uj, Si), as given by Eqs. 2 5 To obey causality, the Green's function corresponding 



where the LDOS p{uj, Si) - 

to ptyp(w) must be specified by analytical continuation, whTcE'is performed by the Hilbert transform 



Gtyp(w) = 



_^ UJ-UJ' 



(7) 



Finally, we close the self-consistency loop by setting the Green's functions of the effective medium be equal to that 
corresponding to the local order parameter, so that 

Ge,„(c^) = Go{oj - - Gtyp(c.). (8) 

It is important to emphasize that our procedure defined by Eqs. [2]|8]is not specific to the problem at hand. The 
same strategy can be used in any theory characterized by a local self-energy. The only requirement specific to our 
problem is the definition of the TDOS as a local order parameter given by Eqj6] . If we choose the algebraic instead 
of the geometric average of the LDOS, our theory would reduce to the well-known coherent potential approximation 
(CPA) |79j, which produces excellent results for the ADOS for any value of disorder, but finds no Anderson transition. 
Thus TMT is a theory having a character very similar to CPA, with a small but crucial difference - the choice of the 
correct order parameter for Anderson localization. 
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Figure 8: Phase diagram for the "semicircular" model. The trajectories of the mobility edge (full line) and the CPA band edge 
(dashed line) are shown as a function the disorder strength W. 



In our formulation, as in DMFT, all the information about the electronic band structure is contained in the choice 
of the bare DOS pa{i-o)- It is not difficult to solve Eqs. |2][8 numerically, which can be efficiently done using Fast 
Fourier transform methods [36_. We have done so for several model of bare densities of states, and find that most 
of our qualitative conclusions do not depend on the specific choice of band structure. We illustrate these findings 
using a simple "semicircular" model for the bare DOS given by po{uj) = — (2w)^, for which Ao{uj) = Go(tLi)/16 

|36| . Here and in the rest of this paper all the energies are expressed in units of the bandwidth, and the random site 
energies Si are uniformly distributed over the interval [—W/2, W/2]. The evolution of the TDOS as a function of W 
is shown in Fig. [TOj The TDOS is found to decrease and eventually vanish even at the band center a,t W 1.36. For 
W < Wc, the part of the spectrum where TDOS remains finite corresponds to the region of extended states (mobile 
electrons), and is found to shrink with disorder, indicating that the band tails begin to localize. The resulting phase 
diagram is presented in Fig. [S] showing the trajectories of the mobility edge (as given by the frequency where the 
TDOS vanishes for a given W, and the band edge where the ADOS as calculated by CPA vanishes. 



B. Critical behavior 



Further insight in the critical behavior is obtained by noting that near W — Wc it proves possible to analytically 
solve Eqs. |2][8) Here we discuss the the critical exponent of the Anderson metal-insulator transition within the TMT 
model. We will demonstrate that the critical exponent /3 with which the order parameter A" vanishes at the transition 
is, in contradiction to the general expectations [7T| . non-universal in this model. 
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1. Critical behavior in the middle of the hand u = Q 



To start with, let us concentrate at the band center (cj = 0), and expand Eqs. 2|8 in powers of the order parameter 
A". In the Hmit oi u — Q self-consistency equations quantities A, G, and E become purely imaginary, and near the 
critical disorder typical Green's function can be expanded in powers of the parameter A": 



G(w,e,) = iA" 



A" 



(w - £, - A')2 + A" 



(9) 



iA" exp 



j deP^Si) \og[e'i + A 



typ 

112} 



iA"f{A")^iA"{a-bA") 



where 



a = /(O) = exp 



-2 J deP{e) \og\e 
.=p[-2/<,.P(., 

-a j deP{e)2TT5[e) = -27raP(0) 



(10) 



-2A" 



e2 + A 



(11) 



and after trivial algebraic operations our self-consistency equations [2][8] reduce to a single equation for the order 
parameter A" 



A" Z'^* 
A" = ^(a-6A") J ^^po{e)e^de. 



(12) 



Equation 12 shows that near the transition along = direction our order parameter A" vanishes linearly (critical 



exponent (3 — 1) independently of the choice of bare lattice DOS pq. In specific case of semicircular bare DOS, where 

A" = aA" - 6A"2. (13) 

the transition where A" vanishes is found at a — 1, giving W — Wc = e/2 = 1.3591, consistent with our numerical 
solution. Near the transition, to leading order 



PtypiW) 



A" 



(14) 



2. Critical behavior near the band edge ui = Uc 

In order to analytically examine scaling the critical behavior at finite w, for simplicity we focus on a semi-circular 
bare DOS, where self-consistency Eqs[2](8]are greatly simplified 



A' + iA" =i> A" = -Trptyp = ImG{uj, Si) 

-A" 



A"(w) = -exp|y"de,P(e,)ln 
A'(c.) = -H[A"{u:)], 



(w - - A') + A"2 



however as in previous section, we expect the critical exponent to be the same for any bare DOS. 

To find the general critical behavior near the mobility edge, we need to expand Eq. 16 in powers of A'' 



A" = A" exp |- y de.P(eO In [(c. - - A')' + A"^] | = A'7(A"), 



(15) 
(16) 
(17) 



(18) 



which cannot be done explicitly, since A' and A" are related via Hilbert transform, which depends on the entire 
function A"(w), and not only on its form near w — lOc- Nevertheless the quantity w — A' {lo) assumes a well-defined 
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W-dependent value at the mobility edge utc = Ldc{W), making it possible for us to determine a range of values a 
critical exponent /3 may take. 



After expanding /(A") = 1 defined by Eq. 18 



/(A") = a-6A" + 0(A"2) 

a = /(0)=exp|-2 J deP(e) In |w - e - A'(w) 



b ^ 



df 



dA" 



= a lim 

A"^o 



A"=0 



deP{e)- 



2A" 



(cj-£- A')2 + A 
= a J deP{e)2TTS{u} - e - A') = 27raP(w - A'), 
we find that to the leading order A" has the following cj dependence {6a = 1 — a) 



(19) 



A" = 



1 



27rP{u} ~ A') 
1 



- 1 



27rP(c^, - A'(c^,)) 
The functional form of (5a (w) is readily found 



6a(u)) cx du)^ . 



a = exp|-2 J deP{e)\n\uj - e - A'{uj) 



Sa{u;) = 2 / rf£P(£)_-_^^(5c. - 5A'(c.)), 



and combining Eqs. 20|22 



we arrive to 



A" = A'^iSuj - SA') 



A" — 



1 



ttP{uj, - A'{uj,)) 



de 



Pie) 



- e - A'(wc) 



(20) 

(21) 
(22) 

(23) 



Note that Suj is negative, since in the range of interest uj < ujc- In Eq 



23 



Ids 



Pie) 

Mw — V is the Hilbert transform of 



P{£), which is positive for uj^ — A' > (right band edge), and it is negative for the left one, where Suj > 0. 

A A" 

A' 

6A'>0 



P=l/4 


5A'<0 







(3=3/4 




Figure 9: Requirement on SA' to be positive definite forces the value of critical exponent /3 to be larger than 1/2 



The lower bound on critical exponent /3 is 0, to insure that A" is convergent and vanishing at a; = coc- Now, if we 
were to assume that the leading contribution to A" comes from Suj (and 6 A' can be neglected), the conclusion would 
be that A" oc Slo, and the critical exponent /3 — 1, just like in w = case. However, this value of /3 is unphysical, 
since the Kramers-Kroning predicts A' to be logarithmically divergent (A' ^ Suj) when A" c>c Soj. This is in direct 
contradiction with our initial statement of Suj being a leading contribution in Eq. 23 (|5A'| <^ \S\), and we conclude 
that 6 A' (X Suj^ is the leading contribution 



A" 



-A'^SA' cx -Sujl^ 



(24) 
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with /3 G (0, 1). This A" being a negative definite quantity imposes a constraint SA' > 0, which is only satisfied for 
/? > 1/2 (see Fig. |9] thus our critical exponent can vary in the range /3 G (1/2, 1). 

Although general arguments for second-order phase transitions |7T] predict universality of exponent /3, we find the 
exponent is non-universal, which is not uncommon in some special cases of mean field theories[80'. It plausible that 
this critical exponent anomaly can be remedied if the MFT is extended to incorporate long range fluctuations effects 
beyond mean-field theory, but this remains an open problem for future work. 
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(0 

Figure 10: Typical density of states for for the SC model, for disorder values W = 0, 0.25, 0.5, 0.75, 1, 1.25, 1.275, 1.3, 1.325, 
1.35. The entire band localizes at W = Wc ^ e/2 ^ 1.359. 




3. Scaling behavior near the critical disorder W — Wc 



The complete analytical solution for TDOS is difficult to obtain for arbitrary cj and W. Still, the approach discussed 
can be extended to find a full frequency-dependent solution ptyp('^j close to the critical value of 



IVBl 



m section 

disorder W = Wc and which assumes a simple scaling form 



(25) 




Figure 11: Left: unsealed typical DOS for various disorder displays parabolic behavior near the MIT. Right: scaling behavior 
near the critical disorder. The range of disorders where parabolic behavior is observed is, in fact, quite broad - W £ (l,Wc), 
VFc = e/2. 
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log(6W) 

Figure 12: Scaling parameters po and cjq as functions of the distance to the transition SW — Wc — W. Numerically obtained 
values (o and •) are in good agreement with analytical predictions (full line). 



distribution of disorder 



A" = -exp 



deP(e) log 



[(w - e - A')2 + A"2] 



after averaging over disorder takes the following form 



2W = a_log(A"2 + a^) + a+log(A"2 + a^ 



-2A" 



arctan 



(^)+arctan(^ 



(28) 



where 



W 

a± = -±(A'-c.). 



(29) 
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Exact expression for the real part of the cavity field A' is obtained by performing a Hilbert transformation of ansatz 



A" = -Trpo 1 - 



a;5 

A' = -i/[A"]=Pof2^-(l-'4)log 



(30) 



U} — UJQ 



Expanding Eqs. |28|30| to the second order in small u results in a system of equations: 

27rpo 



W 
27rpo 



f W \ 1 



4po 



4 

- 1 



arctan = — ^777= — -, 

\27rpoJ 2(i^+^2^2)' 



which can be solved for scaling parameters used the in original ansatz, Eq. |26[ 

4{Wc - W) 



Po 



(31) 
(32) 



Numerical test of TMT 



In order to gauge the quantitative accuracy of out theory, we have carried out first-principles numerical calculations 
for a three dimensional cubic lattice with random site energies, using exact Green functions for an open finite sample 
attached to two semi-infinite clean leads [;41|. We computed both the average and the typical DOS at the band center 
as a function of disorder, for cubes of sizes L =4, 6, 8, 10, 12, and 16, and averages over 1000 sample realizations, in 
order to obtain reliable data by standard finite size scaling procedures. The TMT and CPA equations for the same 
model were also solved by using the appropriate bare DOS (as expressed in terms of elliptic integrals), and the results 
are presented in Fig. [T3| 




Figure 13: Typical and average DOS for a tliree dimensional cubic lattice at the band center (cj — 0). Results from first- 
principle numerical calculations (circles) are compared to the predictions of TMT (for TDOS - full line ) and CPA (for ADOS 
- dashed line). 
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We find remarkable agreement between these numerical data [H] and the self-consistent CPA calculations for the 
ADOS, but also a surprisingly good agreement between the numerical data and the TMT predictions for the TDOS 
order parameter. For a cubic lattice, the exact value is Wc ~ 1.375 [HI], whereas TMT predicts a 20% smaller value 
Wc ~ 1.1. The most significant discrepancies are found in the critical region, since TMT predicts the order parameter 
exponent /3 = 1, whereas the exact value is believed to be /3 « 1.5, consistent with our numerical data. Argument 
based on the multi-fractal scaling analysis |77}l78|. together with numerical calculations[81_ of the multi-fractal spectra 
of wavefunction amplitudes have suggested that in three dimensions, the TDOS order parameter exponent f3 should 
be equal to the conductivity exponent fi 1.5. The result /3 = fi = 1 + 0{e) is also found within the 2 + e approach 
|77L 175] . Nevertheless, we conclude that TMT is as accurate as one can expect from a simple mean-field formulation. 



D. Transport properties 

Most previous conventional transport theories, while providing a wonderful description of good metals, fail to 
describe the transport properties of highly disordered materials. 

In most metals, the temperature coefficient of resistivity (TCR) a is positive, because phonon scattering decreases 
the electronic mean free path as the temperature is raised. The sign of TCR 

d\npres{T) 

a = (33) 

can be deduced from Matthiessen's rule which asserts that the total resistivity in the presence of two or more scattering 
mechanisms is equal to the sum of the resistivities that would result if each mechanism were the only one operating, 
for example: 

Pres= pill +411 (34) 

Matthiessen's rule, as stated in Eq [34] follows from the Boltzmann equation with the assumption of a wave- vector- 
independent relaxation time for each scattering mechanism, so if po is the resistivity of a disordered metal at zero 
temperature and pph(T) is the resistivity of the ordered material due to electron-phonon scattering, then the total 
resistivity at finite temperature is 

p{T) = Pa + pphiT) > po (35) 

predicting that the TCR is positive {a > 0) 

Mooij |H1| in 1973 have pointed out that there exist many highly disordered metals which are poor conductors 
and have p{T) < po and negative TCRs, which clearly violate Matthiessen's rule. In fact in these materials the 
Boltzmann-equation formalism itself is breaking down. Apparently what is happening is that, because of the strong 
disorder and resulting multiple correlated scattering, the Boltzmann hypothesis of independent scattering events fails. 
The simple picture of temperature fluctuations impeding transport of electrons (positive TCR) is now replaced with 
the temperature fluctuations releasing the localized electrons and increasing the conductivity (negative TCR). When 
the transport properties are addressed within the TMT, the interplay of several localization mechanisms is considered, 
which is capable of producing the negative TCRs observed in highly disordered materials. 

We start addressing the transport properties of our system within the TMT by pointing out that the escape rate 
from a given site can be rigorously deflned in terms of the cavity fleld (see Eq. |2| , and using our solution of the TMT 
equations, we flnd t~J = — ImA(O) ^ Ptyp ~ {Wc — W). To calculate the conductivity within our local approach, we 
follow a strategy introduced by Girvin and Jonson (GJ) [S^, who pointed out that close to the localization transition, 
the conductivity can be expressed as cr = Aai2, where A is a vertex correction that represents hops to site outside of 
the initial pair i and j, and 012 is a two-site contribution to the conductivity, that can be expressed as 

ai2 =< A12A21 - A11A22 >, (36) 

where Aij — — ImGy is the spectral function corresponding to the nearest neighbor two-site cluster, < • • • > represents 
the average over disorder. 

We examine the temperature dependence of the conductivity as a function of W . Physically, the most important 
effect of finite temperatures is to introduce finite inelastic scattering due to interaction effects. At weak disorder, 
such inelastic scattering increases the resistance at higher temperatures, but in the localized phase it produces the 
opposite effect, since it suppresses interference processes and localization. To mimic these inelastic effects within our 
noninteracting calculation, we introduce by hand an additional scattering term in our self-energy, viz. S — )• E — 177 or 
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Figure 14: Conductivity as a function of the inelastic scattering rate rj for for the SC model at the band center and W = 0, 
0.125, 0.25, 0.5, 0.75, 1, 1.25, 1.36, 1.5, 1.75, 2. The "separatrix" (cr — a* independent of rj, i. e. temperature) is found at 
W = W* « 1 (dashed Hne). The critical conductivity (Tc(i]) ~ rj^^'^ corresponds to W = Wc = 1.36 (heavy full line). 



it can be treated as the imaginary part of w — w + firy. The parameter r] measures the inelastic scattering rate, and 
is generally expected to be a monotonically increasing function of temperature. 
The relevant ry-dependent Green's functions Gij 



Lu + irj - 



[uj + iri 



{ijj + irj — Ei — A)(cj + irj — Sj — t^) — 
t 



A){uj + iT]- 



t2) - i2 ' 



(37) 



reduce expression |36| to an integrable form 

(A" - ?7)2 
ai2 = 4^ — ...„ X 



dx 



-w/2 x2 + (A"-,7); 



arctan 



-X + y{x^ + (A" - 7?)2) 



(A" - r/)(a;2 + (A" - ry)2 + 1) 



y=L0-A'+W/2 
y=ui-A'-W/2 



(38) 



that can be solved numerically as a function of temperature rj and disorder W . 

We have computed ai2 by examining two sites embedded in the effective medium defined by TMT (Axmt)i thus 
allowing for localization effects. The vertex function A remains finite at the localization transition [84|, and thus can 
be computed within. We have used the CPA approach to evaluate the vertex function as A = Ocpal where Ccpa 
is the CPA conductivity calculated using approach described by Elliot [79] 

f-S/2 



cr(w) cx 



B/2 

Po(e) 



depQ{e)Im[G{ijj, e)] 



52 



3/2 



(39) 
(40) 



and a^2° is the two-site correlation function embedded in the CPA effective medium(AcpA)- Since TMT reduces to 
CPA for weak disorder, our results reduce to the correct value at <C where the conductivity reduces to the 



17 



Drude-Boltzmann form. The resulting critical behavior of the T = conductivity follows that of the order parameter, 
a ptyp {Wc — W), giving the conductivity exponent /i equal to the order parameter exponent /3, consistent with 
what is expected. 



The resulting dependence of the conductivity as a function of 77 and W is presented in Fig. 14 As r] (i. e. 
temperature) is reduced, we find that the conductivity curves "fan out", as seen in many experiment close to the MIT 
l68] . Note the emergence of a "separatrix" |7l [68] where the conductivity is temperature independent, which is found 



for W ~ 1, corresponding to kpi ^ 2, consistent with some experiments [7]. At the MIT, Uc{t]) ~ PiYviv) ~ ?7^^^. 



V. MOTT-ANDERSON TRANSITIONS 



A. Two-fluid picture of Mott 



A first glimpse of the basic effect of disorder on the Mott transition at half filling was outlined already by Mott, [T] 
who pointed out that important differences can be seen even in the strongly localized (atomic) limit. 



w -^^-------^ u 
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Figure 15: Energy level occupation in the strongly localized (atomic) limit. In a Mott-Anderson insulator (center), the disorder 
strength W is larger then the Coulomb repulsion U , and a two-fluid behavior emerges. Here, a fraction of localized states are 
doubly occupied or empty as in an Anderson insulator. Coexisting with those, other states remain singly occupied forming 
local magnetic moments, as in a Mott insulator. 



For weak to moderate disorder W < U , the Mott insulator survives, and each localized orbital is singly occupied 
by an electron, forming a spin 1/2 magnetic moment. For stronger disorder {W > U) the situation is different. Now, 
a fraction of electronic states are either doubly occupied or empty, as in an Anderson insulator. The Mott gap is 
now closed, although a finite fraction of the electrons still remain as localized magnetic moments. Such a state can 
be described |39l HO] as an inhomogeneous mixture of a Mott and an Anderson insulator. A very similar "two-fluid 
model" - of coexisting local magnetic moments and conduction electrons - was proposed |85) 186] some time ago 
on experimental grounds, as a model for doped semiconductors. Some theoretical basis of such behavior has been 
discussed |37fHni [7H [57H55] . but the corresponding critical behavior remains a puzzle. 

This physical picture of Mott is very transparent and intuitive. But how is this strongly localized (atomic) limit 
approached when one crosses the metal-insulator transition from the metallic side? To address this question one needs 
a more detailed theory for the metal-insulator transition region, which was not available when the questions posed by 
Mott and Anderson were first put forward. 



B. Mott or Anderson... or both? 



Which of the two mechanisms dominates criticality in a given material? This is the question often asked when 
interpreting experiments, but a convincing answer is seldom given. To answer it precisely, one must define the 
appropriate criteria - order parameters - characterizing each of the two routes. The conceptually simplest theoretical 
framework that introduces such order parameters is given by TMT-DMFT - which we introduced in the preceding 
section, and discussed in detail in the noninteracting limit.. As in conventional DMFT, its self-consistent procedure 
formally sums-up all possible Feynman diagrams providing local contributions to the electronic self-energy |36) . 
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When the procedure is applied to systems with both interactions and disorder systems, the self energy is still local, 
but may display strong-site-to-site fluctuations. Its low-energy form 

Si(a;„) = (1 - Z,^^) oj„ + - + Ai, 

defines local Fermi liquid parameters |39l 175] : the local quasi-particle (QP) weight Zi, and the renormalized disorder 
potential Vi. This theory portrays a picture of a spatially inhomogeneous Fermi liquid, and is able to track its evolution 
as the critical point is approached. 

In this language, Anderson localization, corresponding to the formation of bound electronic states, is identified by the 
emergence of discrete spectra [31] in the local density of states (LDOS). As we have seen above, this corresponds 1391^1] 
to the vanishing of the typical (geometrically averaged) LDOS ptyp — exp < ln(/9i) >. In contrast, Mott localization 
of itinerant electrons into magnetic moments is identified by the vanishing of the local QP weights (Zj — >■ 0). It is 
interesting and important to note that a very similar physical picture was proposed as the key ingredient for "local 
quantum criticalitv"|90|. or "deconfined quantum criticality" [91', at the T = magnetic ordering in certain 
heavy fermion systems. A key feature in these theories is the possibility that Kondo screening is destroyed precisely 
at the quantum critical point. As a result, part of the electrons - those corresponding to tightly bound f-shells of 
rare earth elements - "drop out" from the Fermi surface and turn into localized magnetic moments. For this reason, 
it is argued, any weak-coupling approach must fail in describing the critical behavior. This is the mechanism several 
groups have attributed to the breakdown of the Hertz-Millis theory |93[ l94] of quantum criticality, which at present 
is believed to be incomplete. 

Precisely the same fundamental problem clearly must be addressed for the Mott-Anderson transition. The trans- 
mutation of a fraction of electrons into local magnetic moments again can be viewed as the suppression of Kondo 
screening - clearly a non-perturbative strong correlation effect - that should be central to understanding the 
critical behavior. To properly characterize it, one must keep track of the evolution of the entire distribution P{Zi) of 
local quasi-particle weights - which can be directly obtained from TMT-DMFT approach |41j to the Mott-Anderson 
transition, which we outlined above. The first applications of this new method to correlated systems with disorder 
was carried out in recent studies by VoUhardt and collaborators [951 196| , who numerically obtained the phase diagram 
for the disordered Hubbard model at half-filling, and discussed the influence of Mott-Anderson localization on mag- 
netically ordered phases. However, the qualitative nature of the critical behavior in the Mott-Anderson transition in 
this model has not been examined in these studies. 



C. Slave-boson solution 

In the following we use complementary semi-analytical methods supplemented by Fermi-liquid theorems, in order 
to clarify the precise form of criticality in this model [97] . By making use of scaling properties |98l 199 1 of Anderson 
impurity models close to the MIT, we present a detailed analytic solution for this problem, which emphasizes the 
dependence of the system properties on its particle-hole symmetry. We consider a half-flUed Hubbard model [40 with 
random site energies, as given by the Hamiltonian 

H = -V ^ cl^Cja + ^ Sirita + ni^n.ii. (41) 

<ij>(T i{7 i 

Here, c^^ (cio-) creates (destroys) a conduction electron with spin tr on site i, riic = cl^Cia, V is the hopping amplitude, 
and U is the on-site repulsion. The random on-site energies Si follow a distribution P(e), which is assumed to be 
uniform and have width W. 

TMT-DMFT [41. ,95 maps the lattice problem onto an ensemble of single-impurity problems, corresponding to sites 
with different values of the local energy e^, each being embedded in a typical effective medium which is self-consistently 
calculated. In contrast to standard DMFT [31], TMT-DMFT determines this effective medium by replacing the 
spectrum of the environment ("cavity") for each site by its typical value, which is determined by the process of 
geometric averaging. For a simple semi-circular model density of states, the corresponding bath function is given by 
[UlinS] A(u;) = V^Gtypiu), with Gtypiuj) = JZ, dw ' Ptyp{i^')/{<^—'-^') being the Hilbert transform of the geometrically- 
averaged (typical) local density of states (LDOS) ptyp{u!) = exp{/ deP{e) In p{uj,e)}. Given the bath function A{uj), 
one flrst needs to solve the local impurity models and compute the local spectra p{io,e) = — tt^^ ImG(w,e), and the 
self-consistency loop is then closed by the the geometric averaging procedure. 

To qualitatively understand the nature of the critical behavior, it is useful to concentrate on the low-energy form 
for the local Green's functions, which can be specifled in terms of two Fermi liquid parameters as 
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where Zi is the local quasi-particle (QP) weight and is the renormalized site energy The parameters Zi and 
£i can be obtained using any quantum impurity solver, but to gain analytical insight here we focus on the variational 
calculation provided by the "four-boson" technique (SB4) of Kotliar and Ruckenstein |100] . which is known to be 
quantitatively accurate at T = 0. The approach consists of determining the site-dependent parameters Ci, di and 
by the following equations 

- ^ - ^ A(w„)G'i(w„) = Z^{^l^e^- e^) ei, (43) 
37 1 

-^-pY. AK)G,(c^„) ^Z,{U-ii-£,+ £,) d„ (44) 

^5]G,(c.„) = iz,(l-e2 + d2), (45) 

where Zi = 2{ei + d,i)^[l — (e^ -I- rf^)]/[l — (ef — df)^] in terms of and di and ^ = U/2. We should stress, though, 
that most of our analytical results rely only on Fermi liquid theorems constraining the qualitative behavior at low 
energy, and thus do not suffer from possible limitations of the SB4 method. 

Within this formulation, the metal is identified by nonzero QP weights Zi on all sites and, in addition, a nonzero 
value for both the typical and the average [pavi'^) — J deP{e)p{ijj,e)] LDOS. Mott localization (i.e. local moment 
formation) is signaled by Zi — > [31], while Anderson localization corresponds to 7^ and pav 7^ 0, but ptyp — 
|32L HT]. While Ref. [53] concentrated on ptyp and pav, we find it useful to simultaneously examine the QP weights 
Zi, in order to provide a complete and precise description of the critical behavior. 



D. Phase diagram 
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Figure 16: T = phase diagram for the disordered half filled Hubbard model, obtained from the numerical SB4 solution of 
TMT-DMFT. 



Using our SB4 method, the TMT-DMFT equations can be numerically solved to very high accuracy, allowing 
very precise characterization of the critical behavior. In presenting all numerical results we use units such that the 
bandwidth B = AV = 1. Fig. 16 shows the resulting T = phase diagram at half filling, which generally agrees 
with that of Ref. [95]. By concentrating first on the critical behavior of the QP weights Zi, we are able to clearly and 
precisely distinguish the metal from the insulator. We find that at least some of the Zi vanish all along the phase 
boundary. By taking a closer look, however, we can distinguish two types of critical behavior, as follows. 
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1. Mott-Anderson vs. Mott-like transition 



For sufficiently strong disorder {W > U), the Mott-Anderson transition proves qualitatively different than the 
clean Mott transition, as seen by examining the critical behavior of the QP weights Zi — Z{ei). Here — >■ only 
for < \ei\ < U/2 , indicating that only a fraction of the electrons turn into localized magnetic moments. The 
rest show Zi 1 and undergo Anderson localization (see below). Physically, this regime corresponds to a spatially 
inhomogeneous system, with Mott fluid droplets interlaced with regions containing Anderson-localized quasiparticles. 
In contrast, for weaker disorder {W < U) the transition retains the conventional Mott character. In this regime 
Zi — > on all sites, corresponding to Mott localization of all electrons. We do not discuss the coexistence region found 
in Ref. because we focus on criticality within the metallic phase. We do not find any "crossover" regime such as 
reported in Ref. [95_ the existence of which we believe is inconsistent with the generally sharp distinction between a 
metal and an insulator at T = 0. 



2. Two-fluid behavior at the Mott-Anderson transition 



To get a closer look at the critical behavior of the QP weights Zi — Z{ei), we monitor their behavior near the 
transition. The behavior of these QP weights is essentially controlled by the spectral weight of our self-consistently- 
determined TMT bath, which we find to vanish at the transition. An appropriate parameter to measure the distance 
to the transition is the bandwidth t of the bath spectral function, which is shown in Fig. 17. 
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Figure 17: Frequency dependence of the typical DOS very close to the metal-insulator transition for (a) the Mott-Anderson 
transition {W > U) at U = 1.25 and (b) the Mott-like transition (W < U) at W — 1.0. The insets show how, in both cases, 
the ptyp{u}) bandwidth t — > at the transitions. 



Considering many single-impurity problems, we observe a two- fluid picture, just as in the limit earlier analyzed by 
Mott. [T] Indeed, these results correspond to the same atomic limit discussed by Mott, since, although the hopping 
itself is still finite, the cavity field "seen" by the impurities goes to zero in the current case. 

As in the atomic limit, the sites with le^j < U/2 turn into local moments and have vanishing quasiparticle weight 
Zi 0. The remaining sites show ^ 1, as they are either doubly occupied, which corresponds to those with 
Ei < —U/2, or empty, which is the case for those sites with Si > U/2. Consequently, as the transition is approached, 
the curves Z{ei,t) "diverge" and approach either Z = or Z = 1. These values of Z can thus be identified as two 
stable fixed points for the problem in question, as we discuss below. 
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Figure 18: Quasiparticle weight Z plotted as a function of the distance to the Mott-Anderson transition for different values 
of the local site energy tjW . We present the results only for positive site energies, as a similar behavior holds for negative ones. 

Note that in Fig. [18] we restrict the results to positive energy values, as a similar behavior is observed for negative 
Ei. In this case, there is precisely one value of the site energy Si = e*, for which Z{£* ,t) — Z* . This corresponds 
to the value of £i below which Z "flows" to and above which Z "flows" to 1. In other words, it corresponds to an 



unstable fixed point. Just as in the atomic limit, e* is equal to U /2 {e* /W — 0.3125 in Fig. 18 where U — 1.75 and 
W = 2.8). 
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Figure 19: Quasiparticle weight Z as a function of t/t*{Se) showing that the results for different e can be collapsed onto a single 
scaling function with two branches. The results for different e correspond to different symbols. The inset shows the scaling 
parameter t* as a function of e/W for the upper (squares) and bottom (circles) branches. 
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3. /3-function formulation of scaling 

Our numerical solutions provide evidence that as a function of t the "charge" Z{t) "flows" away from the unstable 
"fixed point" Z* , and towards either stable "fixed points" Z — or Z — 1. The structure of these fiows show power-law 
scaling as the scale t — )■ 0; this suggests that it should be possible to collapse the entire family of curves Z{t, Se) onto 
a single universal scaling function 

Zit,Se) = f[t/t*{Se)l (46) 
where the crossover scale t*{6e) = C^\6£\'^ around the unstable fixed point. Remarkably, we have been able to scale 



the numerical data precisely in this fashion, see Fig. 19 and extract the form of t*{Se). We find that t*(6s) vanishes 
in a power law fashion at 5e = 0, with exponent (f) — 2 and the amplitudes differ by a factor close to two for 
Z^Z*. 

As shown in Fig. 19 the scaling function f{x) where x ~ t/t*{5e) presents two branches: one for < e* and 
other for £i > e* . We found that for x — both branches of f{x) are linear in x, while for x ^ 1 they merge, i.e. 
f{x) — > Z* ± A'^x~^/'^. As can be seen in the first two panels, in the limit t — 0, the curve corresponding to < U /2 
has Z{t) = B^t, while that for Ei > U/2 follows 1 — Z{t) = B^t. These results are for a flat cavity field but, as 
mentioned earlier, we checked that the same exponents are found also for other bath functions, meaning that they 
are independent of the exact form of the cavity field. The power-law behavior and the respective exponents observed 
numerically in the three limits above have also been confirmed by solving the SB equations analytically |99|close to 
the transition (t — 0) . 

In the following, we rationalize these findings by defining an appropriate /^-function which describes all the fixed 
points and the corresponding crossover behavior. Let us assume that 

is an explicit function of Z only, but not of the parameters t or 6e. The desired structure of the fiows would be 
obtained if the /3-function had three zeros: at Z = and Z = 1 with negative slope (stable fixed points) and one at 
Z — Z* with positive slope (unstable fixed point). The general structure of these fiows can thus be described in a 
/3-function language similar to that used in the context of a renormalization group approach; we outline the procedure 
to obtain /3(Z) from the numerical data. 




Figure 20: /3-function obtained as described in the text for the Anderson impurity models close to the Mott- Anderson transition. 
The filled circles indicate the three fixed points found for this problem. The arrows indicate how Z flows to the stable points 
{Z = and Z — 1) and from the unstable one (Z ~ 0.7). 
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The integration of Eq. (47 1 can be written in the form of Eq. (46 I as 



Z = f[t/t*{Zo)] 



where Zo is the initial condition (a function of Se). With x = t/t* as before, Eq. (46) can be rewritten as 

= -xf{x). 



(48) 



(49) 



The numerical data for Z = f{x) as a function of x is presented in Fig. 19 Thus, using Eq. (49 1, the /3-function 
in terms of x{Z) is determined, which can finally be rewritten in terms of Z. Carrying out this procedure, we obtain 
/3{Z) as shown in Fig. 20 In accordance with what was discussed above, we see that fi{Z) has three fixed points, as 
indicated in the figure by filled circles. Z = and Z = 1 are stable, while Z « 0.7 is the unstable fixed point. The 
scaling behavior and the associated /3-function observed here reflect the fact these impurity models have two phases 
(singlet and doublet) when entering the insulator. The two stable fixed points describe these two phases, while the 
unstable fixed point Z* describes the phase transition, which is reached by tuning the site energy. 

Interestingly, the family of curves in Fig. [18] looks similar to those seen in some other examples of quantum critical 
phenomena. In fact, one can say that the crossover scale t plays the role of the reduced temperature, and the reduced 
site energy Se ~ (si — £*) /s* that of the control parameter of the quantum critical point. As the site energy is 
tuned at t = 0, the impurity model undergoes a phase transition from a singlet to a doublet ground state. Quantum 
fluctuations associated with the metallic host introduce a cutoff and round this phase transition, which becomes sharp 
only in the i — >■ limit. 



E. Wavefunction localization 



To more precisely characterize the critical behavior we now turn our attention to the spatial fluctuations of the 
quasiparticle wavef unctions, we compare the behavior of the typical {ptyp) and the average {pav) LDOS. The approach 
to the Mott-Anderson transition {W > U) is illustrated by increasing disorder W for fixed U = 1.25 (Fig. 21 - top 
panels). Only those states within a narrow energy range (w < see also Fig. 17) around the band center (the Fermi 
energy) remain spatially delocalized {ptyp ~ Pav), due to strong disorder screening |3H[99] within the Mott fiuid (sites 
showing Zi ^ at the transition). The electronic states away from the band center (i.e. in the band tails) quickly get 
Anderson-localized, displaying large spatial fluctuations of the wavefunction amplitudes |41| and having ptyp ^ Pav 
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Figure 21: Frequency dependence of ptyp (full line) and pav (dashed line) in the critical region. Results in top panels illustrate 
the approach to the Mott-Anderson transition {W > U) ai U — 1.25; the bottom panels correspond to the Mott-like transition 
{W < U) at W = 1.0. For the Mott-Anderson transition, only a narrow band of delocalized states remain near the Fermi energy, 
corresponding to ptyp 7^ 0. In contrast, most electronic states remain delocalized ptyp ~ pav near the Mott-like transition. 
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The spectral weight of the delocaHzed states (states in the range u> < t) decreases with disorder and vanishes at the 
transition, indicating the Mott localization of this fraction of electrons. At this critical point, the crossover scale t also 
vanishes. In contrast, the height ptyp(0) remains finite at the transition, albeit at a reduced VK-dependent value, as 
compared to the clean limit. More precise evolution of ptyp{0) is shown in Fig. 22a, demonstrating its critical jump. 

Behavior at the Mott-like transition {W < U) is dramatically different (Fig. 21 - bottom panel). Here ptyp « Pav 
over the entire QP band, indicating the absence of Anderson localization. It proves essentially identical as that 
established for the disordered Hubbard model within standard DMFT |31,, reflecting strong correlation-enhanced 
screening of disorder [31, 99,, where both Pau(w = 0) and ptyp{uJ = 0) approach the bare {W = 0) value (see also 
Fig. 22b). Similar results were found in Ref. [53], but an explanation was not provided. 

The corresponding pinning |3H [99] for p{uj = 0, e) is shown in the insets of Fig. 3, both for the Mott- Anderson and 
the Mott-like transition. In the Mott- Anderson case, this mechanism applies only within the Mott fluid (|£| < U/2), 
while within the Anderson fluid (|e| > U/2) it assumes smaller values, explaining the reduction of ptyp{0) in this case. 



F. Analytical solution 

Within our SB4 approach, the TMT-DMFT order-parameter function ptyp{uj) satisfles the following self-consistency 
condition 

ptypiuj)=exp J deP{e){\n[V^Z^{e)ptyp{uj)] 
- ln[(cj - e(£) - V^Zie) Re Gtyp{uj) f 

+{nV'Z{e)ptyp{u:)f]}. (50) 

While the solution of this equation is in general difficult, it simplifles in the critical region, where the QP parameter 
functions Z{e) and e{e) assume scaling forms which we carefully studied in previous work [13]. This simplification 
allows, in principle, to obtain a closed solution for all quantities. In particular, the crossover scale t, which defines 
the ptyp{ijj) mobility edge (see Fig. 17 and Ref. [99 ), is determined by setting ptyp{uj — t) = 0. 

Using this approach we obtain that, in the case of Mott-like transition {W < U), the critical behavior of all 
quantities reduces to that found in standard DMFT [31], including t Uc{W) — U (in agreement with the numerical 
results of Fig. 17b), perfect screening of site randomness |3U WM> ^^'^ the approach of Pav{^ = 0) and ptyp{i^ = 0) 
to the clean value. The precise form of the critical behavior for the crossover scale t is more complicated for the 
Mott-Anderson transition (W > U) (as confirmed by our numerical results in Fig. 17a), and this will not be discussed 
here. 

Instead, we focus on elucidating the origin of the puzzling behavior of p^ — ptyp{i-^ = 0), which is known [41 to 
vanish linearly pc ~ {Wc — W) for [/ = 0, but which we numerically find to display a jump (i.e. a finite value) at 
criticality, as soon as interactions are turned on. For lu = our self-consistency condition reduces (for our model 
ReGtyp(O) = by particle-hole symmetry) to 

deP (e) In ..^r^ n 9 = 0, (51) 

which further simplifies as we approach the critical point. Here, the QP parameters Z{e) — > and e{e) ~ Z'^is) ^ 
Z{e) for the Mott fluid (|e| < U/2), while Z{e) — > 1 and |£(e)| — ^ |e - U/2\ for the Anderson fluid (|£| > U/2), and 
we can write 

fU/2 



0= / d£P(£)ln 

(W-U)/2 



2 



deP (£) ln[(£/l/)' + {nVp,)% (52) 





This expression becomes even simpler in the U «W limit, giving 

^ln-i-+a-6^p, + 0[p2] =0, (53) 

where a{W,U) = (1 - U /W){1 - ln[(W - U)/2V]} and b = This result reproduces the known result [H] 

Pc ~ (Wc — W) at U = 0, but dramatically different behavior is found as soon as J7 > 0. Here, a non-analytic 



25 



O 



o 
ol 



1.0 
0.5 
0.0 
1.2 
1.1 
1.0 



-0.8 0.0 0.8 1. 
- ' I |_ 1^ ^ I ' 



(a) 



0.5 
G— 



1.0 
"T" 
-a- 



W 
1.5 

"r 



2.0 



2.5 



3.0 



^-aa^ 



o 
1^ 



G O typical 

□ — □ average 



h 
P 
1 

-e- 



O— — O typical 
_ □ — □ average 



h (b) er 



1.3 



^ I 
6 -0.3 0.0 0.3| 

' I ' i ' I ' 



1.251- / 
_ / 



\ 



0.5 



1.0 



1.5 



2.0 



U 



Figure 22: Typical and average values of p(0) as the metal-insulator transition is approached for (a) U — 1.25 and (b) W = 1.0. 
The insets show p(0) as a function of e for (a) W = 2.5, 2.6, 2.7 and 2.83 (from the black curve to the blue one) and (b) 
U = 1.5, 1.6, 1.7 and 1.86. 



(singular) contribution emerges from the Mott fluid (|e| < U/2), which assures that pc must remain finite at the 
critical point, consistent with our numerical results (see Fig. 22). Note that the second term in Eq. (52 1, coming from 
the Anderson fluid {\e\ > U/2), vanishes in the case of a Mott-like transition {U > W), and our result reproduces the 
standard condition npcV = 1 |31j . which corresponds to the clean limit. 

A further glimpse on how the condition wpcV = 1 is gradually violated as we cross on the Mott-Anderson side is 
provided by solving Eq. (52) for U limit, giving 



1 



24 \V 



U 
W 



(54) 



again consistent with our numerical solution. 



But what is the physical origin of the jump in pcl To see it, note that the singular form of the first term in Eq. ( 52 1 



comes from the Kondo pinning j31J e{e) ~ Z^is) ^ Z{e) within the Mott fiuid. This behavior refiects the particle-hole 
symmetry of our (geometrically averaged) ptyp{i^ = 0) bath function, which neglects site-to-site cavity fiuctuations 
present, for example, in more accurate statDMFT theories [291 15U1 [5^ I5U1 [751 [75] . Indeed, in absence of particle-hole 
symmetry, one expects ||31J e^e) ~ Z{£), and the resulting e-dependence should cut-off the log singularity responsible 
for the jump in pc- This observation provides a direct path to further refine the TMT-DMFT approach, reconciling 
the present results with previous statDMFT findings |29| [30| [39l [60j [75| [76] . As a next step, one should apply the 
TMT ideas to appropriately chosen effective models [62] , in order to eliminate those features reflecting the unrealistic 
paxticle-hole symmetry built in the current theory. We emphasize that the two-fluid picture is a consequence of only 
a fraction of the sites showing Z — >■ and is not dependent on either particle-hole symmetry or the consequent jump 
in the DOS. This interesting research direction is just one of many possible future applications of our TMT-DMFT 
formalism. 
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VI. CONCLUSIONS AND OUTLOOK 

This article described the conceptually simplest theoretical approach which is able to capture the interplay of strong 
correlation effects - the Mott physics - and the disorder effects associated with Anderson localization, ft demonstrated 
that one can identify the signatures of both of these basic mechanisms for localization by introducing appropriate local 
order parameters, which are then self-consistently calculated within the proposed Typical-Medium Theory. We 
showed that key insight can be obtained by focusing on the evolution of the local quasiparticle weights as a second 
order parameter describing tendency to Mott localization, in addition to the Anderson-like TMT order parameter 
Ptyp- Our main finding is that, for sufficiently strong disorder, the physical mechanism behind the Mott-Anderson 
transition is the formation of two fluids, a behavior that is surprisingly reminiscent of the phenomenology proposed 
for doped semiconductors HOlJ. Here, only a fraction of the electrons (sites) undergo Mott localization; the rest 
can be described as Anderson-localized quasiparticles. Physically, it describes spatially inhomogeneous situations, 
where the Fermi liquid quasiparticles are destroyed only in certain regions - the Mott droplets - but remain coherent 
elsewhere. Thus, in our picture the Mott-Anderson transition can be seen as reminiscent of the "orbitally selective" 
Mott localization |102[ 1103] . To be more precise, here we have a "site selective" Mott transition, since it emerges in 
a spatially resolved fashion. Understanding the details of such "site selective" Mott transitions should be viewed as 
an indispensable first step in solving the long-standing problem of metal- insulator transitions in disordered correlated 
systems. 

Acknowledgements 

The author thanks Elihu Abrahams, Carol Aguiar, Eric Andrade, Gabi Kotliar, Eduardo Miranda, Andrei Pastor, 
and Darko Tanaskovic for many years of exciting and fruitful collaboration. This work was supported by the NSF 
grant DMR-0542026. 



[1] N. F. Mott, Metal- Insulator Transition. (Taylor & Francis, London, 1990). 

[2] Y. Kohsaka, C. Taylor, K. Fujita, A. Schmidt, C. Lupien, T. Hanaguri, M. Azuma, M. Takano, H. Eisaki, H. Takagi, 
S. Uchida, and J. C. Davis, An intrinsic bond-centered electronic glass with unidirectional domains in underdoped cuprates, 
Science. 315, 1380, (2007). 

[3] W. D. Wise, K. Chatterjee, M. C. Boyer, T. Kondo, T. Takeuchi, H. Ikuta, Z. Xu, J. Wen, G. D. Gu, Y. Wang, and 
E. W. Hudson. Imaging nanoscale Fermi surface variations in an inhomogeneous superconductor, (2008). 

[4] A. N. Pasupathy, A. Pushp, K. K. Gomes, C. V. Parker, J. Wen, Z. Xu, G. Gu, S. Ono, Y. Ando, and A. Yazdani, 
Electronic origin of the inhomogeneous pairing interaction in the high-tc superconductor Bi2Sr2CaCu208+d, Science. 
320, 196, (2008). 

[5] E. Dagotto, Complexity in strongly correlated electronic systems. Science. 309, 257, (2005). 

[6] E. Miranda and V. Dobrosavljevic, Disorder-driven non-Fermi liquid behavior of correlated electrons. Reports on Progress 

in Ptiysics. 68, 2337, (2005). 
[7] P. A. Lee and T. V. Ramakrishnan, Disordered eletronic systems. Rev. Mod. Pfiys. 57(2), 287 (4, 1985). 
[8] E. Abrahams, P. W. Anderson, D. C. Licciardello, and T. V. Ramakrishnan, Scaling theory of localization: Absence of 

quantum diffusion in two dimensions, Phys. Rev. Lett. 42, 673-676, (1979). 
[9] F. Wegner, The mobility edge problem: continuous symmetry and a conjecture, Phys. Rev. B. 35, 783, (1979). 
[10] L. Schaffer and F. Wegner, Disordered system with n orbitals per site: Lagrange formulation, hyperbolic symmetry, and 

goldstone modes, Ptiys. Rev. B. 38, 113, (1980). 
[11] N. Goldenfeld, Lectures on Ptiase Transitions and ttie Renorrnalization Group. (Addison- Wesley, New York, 1992). 
jl2] K. B. Efetov, A. I. Larkin, and D. E. Khmel'nitskii, Interaction between diffusion modes in localization theory, Zh,. Eksp. 
Teor. Fiz. 79, 1120, (1980). 

[13] A. M. Finkel'stein, Influence of Coulomb interaction on the properties of disordered metals, Zti. Eksp. Teor. Fiz. 84, 168, 

(1983). [Sov. Phys. JETP 57, 97 (1983)]. 
114] A. M. Finkel'stein, Zh. Eksp. Teor. Fiz. 86, 367, (1984). ]Sov. Phys. JETP 59, 212 (1983)]. 

]15] C. Castellani, C. D. Castro, P. A. Lee, and M. Ma, Interaction-driven metal-insulator transitions in disordered fermion 

systems, Phys. Rev. B. 30, 527, (1984). 
116] D. Belitz and T. R. Kirkpatrick, The Anderson-Mott transition. Rev. Mod. Phys. 66(2), 261 - 380 (Apr., 1994). 
117] A. Punnoose and A. M. Finkel'stein, Dilute electron gas near the metal-insulator transition: Role of valleys in silicon 

inversion layers, Phys. Rev. Lett. 88, 016802(4), (2002). 
118] C. Castellani, B. G. Kotliar, and P. A. Lee, Fermi-liquid theory of interacting disordered systems and the scaling theory 

of the metal-insulator transition, Phys. Rev. Lett. 56, 1179, (1987). 
119] R. B. Grifliths, Correlations in Ising ferromagnets. I, J. Math. Phys. 8(3), 478 (Mar., 1967). 



27 



R. B. Griffiths, Nonanalytic behavior above the critical point in a random Ising ferromagnet, Phys. Rev. Lett. 23, 17-19 

(1 July, 1969). 

H. V. Lohneysen, A. Rosch, M. Vojta, and P. Wolfle, Fermi-liquid instabilities at magnetic quantum phase transitions. 
Rev. Mod. Phys. 79(3), 1015 (17 Aug., 2007). 

S. Sachdev, Quantum Phase Transitions. (Cambridge University Press, UK, 1999). 

T. Vojta, Rare region effects at classical, quantum and nonequilibrium phase transitions, J. Phys. A: Math. Gen. 39, 
R143-R205 (16 May, 2006). 

D. S. Fisher, Critical behavior of random transverse-field Ising spin chains, Phys. Rev. B. 51(10), 6411-6461 (1 Mar., 

1995). 

T. Vojta, Disorder-induced rounding of certain quantum phase transitions, Phys. Rev. Lett. 90, 107202 (12 Mar., 2003). 
J. A. Hoyos and T. Vojta, Theory of smeared quantum phase transitions, Phys. Rev. Lett. 100(24), 240601 (17 June, 
2008). 

V. Dobrosavljcvic and E. Miranda, Absence of conventional quantum phase transitions in itinerant systems with disorder, 
Phys. Rev. Lett. 94, 187203 (11 May, 2005). 

M. J. Case and V. Dobrosavljevic, Quantum critical behavior of the cluster glass phase, Phys. Rev. Lett. 99(14), 147204 
(3 Oct., 2007). 

E. Andrade, E. Miranda, and V. Dobrosavljevic, Energy-resolved spatial inhomogeneity of disordered mott systems, 
Physica B: Condensed Matter. 404(19), 3167 - 3171, (2009). 

E. C. Andrade, E. Miranda, and V. Dobrosavljevic, Electronic Griffiths phase of the d=2 Mott transition. Physical Review 
Letters. 102, 206403, (2009). 

D. Tanaskovic, V. Dobrosavljcvic, E. Abrahams, and G. Kotliar, Disorder screening in strongly correlated systems, Phys. 
Rev. Lett. 91, 066603, (2003). 

P.W.Anderson, Absence of diffusion in certain random lattices, Phys. Rev. 109, 1492-1505, (1958). 

N. E. Hussey, K. Takenaka, and H. Takagi, Universality of the Mott-Ioff-Regel limit in metals. Philosophical Magazine. 

84, 2847, (2004). 

M. M. Radonjic, D. Tanaskovic, V. Dobrosavljevic, and K. Haule, Influence of disorder on incoherent transport near the 
mott transition, Phys. Rev. B. 81(7), 075118 (Feb, 2010). 

J. G. Analytis, A. Ardavan, S. J. Blundell, R. L. Owen, E. F. Carman, C. Jeynes, and B. J. Powell, The effect of irradiation- 
induced disorder on the conductivity and critical temperature of the organic superconductor k;-(BEDT-TTF)2Cu(SCN)2, 
Phys. Rev. Lett. 96, 177002, (2006). 

A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, Dynamical mean-field theory of strongly correlated fermion 

systems and the limit of infinite dimensions. Rev. Mod. Phys. 68, 13, (1996). 

V. Dobrosavljcvic, T. R. Kirkpatrick, and B. G. Kotliar, Kondo effect in disordered systems, Phys. Rev. Lett. 69, 1113, 
(1992). 

V. Dobrosavljevic and G. Kotliar, Strong correlations and disorder in d=oo and beyond, Phys. Rev. B. 50, 1430, (1994). 
V. Dobrosavljevic and G. Kotliar, Mean field theory of the Mott-Anderson transition, Phys. Rev. Lett. 78, 3943, (1997). 
V. Dobrosavljevic and G. Kotliar, Dynamical mean-field studies of metal-insulator transitions, Philos. Trans. R. Soc. 
London A. 356, 57, (1998). 

V. Dobrosavljevic, A. Pastor, and B. K. Nikolic, Typical medium theory of Anderson localization: A local order parameter 
approach to strong disorder effects, Europhys. Lett. 62, 76-82, (2003). 

M. C. O. Aguiar, E. Miranda, and V. Dobrosavljevic, Localization effects and inelastic scattering in disordered heavy 
electrons, Phys. Rev. B. 68, 125104, (2003). 

M. C. O. Aguiar, E. Miranda, V. Dobrosavljevic, E. Abrahams, and G. Kotliar, Temperature dependent transport of 

correlated disordered electrons: elastic vs. inelastic scattering, Europhys. Lett. 67, 226, (2004). 

E. Miranda and V. Dobrosavljcvic, Localization effects in disordered Koudo lattices, Physica B. 259-261, 359, (1999). 
P. W. Anderson, Local moments and localized states, Rev. Mod. Phys. 50(2), 191, (1978). 

A. A. Pastor and V. Dobrosavljevic, Melting of the electron glass, Phys. Rev. Lett. 83, 4642, (1999). 

V. Dobrosavljevic, D. Tanaskovic, and A. A. Pastor, Glassy behavior of electrons near metal-insulator transitions, Phys. 

Rev. Lett. 90, 016402, (2003). 

D. Dalidovich and V. Dobrosavljevic, Landau theory of the Fermi-liquid to electron-glass transition, Phys. Rev. B. 66, 
081107, (2002). 

L. Arrachea, D. Dalidovich, V. Dobrosavljevic, and M. J. Rozenberg, Melting transition of an Ising glass driven by a 
magnetic field, Phys. Rev. B. 69(6), 064419 (Feb, 2004). 

S. Pankov and V. Dobrosavljevic, Nonlinear screening theory of the Coulomb glass, Phys. Rev. Lett. 94, 046402, (2005) 
M. MuUer and L. B. loffe, Glass transition and the Coulomb gap in electron glasses, Phys. Rev. Lett. 93, 256403, (2004) 
S. Sachdev and N. Read, Metallic spin glasses, J. Phys. Condens. Matter. 8, 9723, (1996). 

A. M. Sengupta and A. Georges, Non-Fermi-liquid behavior near a T=0 spin-glass transition, Phys. Rev. B. 52, 10295 
(1995). 

H. Westfahl Jr, J. Schmalian, and P. G. Wolynes, Dynamical mean field theory for self-generated quantum glasses, Phys 
Rev. B. 68, 134203, (2003). 

S. Wu, J. Schmalian, G. Kotliar, and P. G. Wolynes, Solution of local-field equations for self-generated glasses, Phys. Rev. 

B. 70(2), 024207 (Jul, 2004). 

V. Dobrosavljcvic, T. R. Kirkpatrick, and G. Kotliar, Phys. Rev. Lett. 69, 1113, (1992). 

E. Miranda, V. Dobrosavljevic, and G. Kotliar, Kondo disorder: a possible route towards non-Fermi liquid behavior, J. 



28 



Phys.: Condens. Matter. 8, 9871, (1996). 
[58] E. Miranda, V. Dobrosavljcvic, and G. Kotliar, Disorder-driven non-Fcrrrii liquid behavior in Kondo alloys, Phys. Rev. 
Lett. 78, 290, (1997). 

[59] E. Miranda, V. Dobrosavljevic, and G. Kotliar, Non-Fermi liquid behavior as a consequence of Kondo disorder, Physica 
B. 230, 569, (1997). 

}60] E. Miranda and V. Dobrosavljevic, Localization-induced Griffiths phase of disordered Anderson lattices, Phys. Rev. Lett. 

86, 264, (2001). 

}61] E. Miranda and V. Dobrosavljevic, Griffiths phase of the Kondo insulator fixed point, J. Magn. Magn. Mat. 226-230, 

110, (2001). 

}62] D. Tanaskovic, E. Miranda, and V. Dobrosavljevic, Effective model of the electronic Griffiths phase, Phys. Rev. B. 70, 
205108, (2004). 

}63] D. Tanaskovic, V. Dobrosavljevic, and E. Miranda, Spin liquid behavior in electronic Griffiths phases, Physical Review 

Letters. 95, 167204, (2005). 
164] B. Shklovskii and A. Efros, Electronic Properties of Doped Semiconductors. (Springer- Verlag, 1984). 
165] R. N. B. M. A. Paalanen, Transport and thermodynamic properties across the metal-insulator transition, Physica B. 169, 

231, (1991). 

166] S. V. Kravchcnko, W. E. Mason, G. E. Bowker, J. E. Furneaux, V. M. Pudalov, and M. D'lorio, Scaling of an anomalous 

metal-insulator transition in a two-dimensional system in silicon at 6 = 0, Phys. Rev. B. 51, 7038, (1995). 
167] D. Popovic, A. B. Fowler, and S. Washburn, Metal-insulator transition in two dimensions: effects of disorder and magnetic 

field, Phys. Rev. Lett. 79, 1543, (1997). 
168] E. Abrahams, S. V. Kravchenko, and M. P. Sarachik, Colloquium: Metallic behavior and related phenomena in two 

dimensions, Rev. Mod. Phys. 73, 251-266, (2001). 
169] D. Simonian, S. V. Kravchenko, and M. P. Sarachik, Reflection symmetry at a b=0 metal-insulator transition in two 

dimensions, Phys. Rev. B. 55(20), R13421-R13423 (May, 1997). 
170] V. Dobrosavljevic, E. Abrahams, E. Miranda, and S. Chakravarty, Scaling theory of two-dimensional metal-insulator 

transitions, Phys. Rev. Lett. 79, 455, (1997). 
171] N. Goldcnfcld, Lectures on phase transitions and the renormalization group. (Addison- Wesley, Reading, 1992). 
172] M. Rozcnberg, G. Kotliar, and H. Kajueter, Transfer of spectral weight in spectroscopies of correlated electron systems, 

Phys. Rev. B. 54, 8542, (1996). 
173] M. Mezard, G. Parisi, and M. A. Virasoro, Spin Glass theory and beyond. (World Scientific, Singapore, 1986). 
174] V. Dobrosavljevic and G. Kotliar, Hubbard models with random hopping in d=oo, Phys. Rev. Lett. 71, 3218, (1993). 
175] V. Dobrosavljcvic and G. Kotliar, Dynamical mean-field studies of metal-insulator transitions, Phil. Trans. R. Soc. Lond. 

A. 356, 1, (1998). 

176] E. Miranda and V. Dobrosavljevic, Disorder-driven non-Fermi liquid behaviour of correlated electrons, Rep. Prog. Phys. 
68, 2337-2408 (22 Aug., 2005). 

177] M.Janssen, Statistics and scaling in disordered mesoscopic electron systems, Phys. Rep. 295, 1, (1998). 

178] A.D.Mirlin, Statistics of energy levels and eigenfunctions in disordered systems, Phys. Rep. 326, 259, (2000). 

179] P. R.J.Elliot, J.A.Krumhansl, The theory and properties of randomly disordered crystals and related physical systems. 
Rev. Mod. Phys. 46, 465, (1974). 

180] H. Sompolinsky and A. Zippelius, Relaxational dynamics of the edwards-andcrsoii iriodel and the rncan-ficld theory of 
spin-glasses, Phys. Rev. B. 25, 6860, (1982). 

181] M. H.Grussbach, Determination of the mobility edge in the anderson model of localization in three dimensions by multi- 
fractal analysis, Phys. Rev. B. 51, 663, (1995). 

182] J. Mooji, Electrical conduction in concentrated disordered transition metal alloys, Phys. Status Solidi A Phys. 17, 521, 
(1973). 

183] M. S.M.Girvin, Dynamical electron-phonon interaction and conductivity in strongly disordered metal alloy, Phys. Rev. 

B. 22, 3583, (1980). 

184] S. M. Girvin and M. Jonson, Dynamical electron-phonon interaction and conductivity in strongly disordered metal alloys, 

Phys. Rev. B. 22(8), 3583-3597 (4, 1980). 
185] J. D. Quirt and J. R. Marko, Phys. Rev. Lett. 26, 318, (1971). 

186] M. A. Paalanen, J. E. Graebner, R. N. Bhatt, and S. Sachdev, Phys. Rev. Lett. 61, 597, (1988). 

187] M. Milovanovic, S. Sachdev, and R. N. Bhatt, Effective-field theory of local-moment formation in disordered metals, Phys. 
Rev. Lett. 63, 82, (1989). 

188] R. N. Bhatt and D. S. Fisher, Absence of spin diffusion in most random lattices, Phys. Rev. Lett. 68, 3072, (1992). 
189] M. Lakner, H. V. Lohneysen, A. Langenfeld, and P. Wolfle, Phys. Rev. B. 50, 17064, (1994). 

190] Q. Si, S. Rabello, K. Ingersent, and J. L. Smith, Locally critical quantum phase transitions in strongly correlated metals. 
Nature. 413, 804, (2001). 

191] T. Senthil, S. Sachdev, and M. Vojta, Fractionalizcd Fermi liquids, Phys. Rev. Lett. 90, 216403, (2003). 

192] T. Senthil, M. Vojta, and S. Sachdev, Weak magnetism and non-Fermi liquids near heavy-fermion critical points, Phys. 

Rev. B. 69, 035111, (2004). 
193] J. A. Hertz, Quantum critical phenomena, Phys. Rev. B. 14, 1165, (1976). 

194] A. J. Millis, Effect of a nonzero temperature on quantum critical points in itinerant fermion systems, Phys. Rev. B. 48 
(10), 7183, (1993). 

195] K. Byczuk, W. Hofstetter, and D. VoUhardt, Mott-Hubbard transition versus Anderson localization in correlated electron 



29 



systems with disorder, Phys. Rev. Lett. 94, 056404, (2005). 
[96] K. Byczuk, W. Hofstetter, and D. VoUhardt, Competition between Anderson localization and antifcrromagnetism in 

correlated lattice fermion systems with disorder. Physical Review Letters. 102, 146403, (2009). 
[97] M. C. O. Aguiar, V. Dobrosavljevic, E. Abrahams, and G. Kotliar, Critical behavior at Mott-Anderson transition: a 

TMT-DMFT perspective, Physical Review Letters. 102, 156402, (2009). 
198] M. C. O. Aguiar, V. Dobrosavljevic, E. Abrahams, and G. Kotliar, Scaling behavior of an Anderson impurity close to the 

Mott-Anderson transition, Physical Review B. 73, 115117, (2006). 
199] M. C. O. Aguiar, V. Dobrosavljevic, E. Abrahams, and G. Kotliar, Disorder screening near the Mott-Anderson transition, 

Physica B. 403, 1417, (2008). 

llOO] G. Kotliar and A. E. Ruckenstein, New Functional Integral Approach to Strongly Correlated Fermi Systems: The 

Gutzwiller Approximation as a Saddle Point, Phys. Rev. Lett. 57(11), 1362, (1986). 
JlOl] M. A. Paalanen, J. E. Graebner, R. N. Bhatt, and S. Sachdev, Thermodynamic Behavior near a Metal-Insulator Transition, 

Phys. Rev. Lett. 61, 597, (1998). 
[102] L. De Leo, M. Civelli, and G. Kotliar, T=0 heavy fermion quantum critical point as an orbital selective Mott transition. 

Physical Review Letters. 101, 256404, (2008). 
1103] C. Pepin, Kondo breakdown as a selective Mott transition in the Anderson lattice, Physical Review Letters. 98, 206401, 

(2007). 



